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VELOCITY STRUCTURE IN THE ORION NEBULA. 
II. EMISSION LINE ATLAS OF PARTIALLY IONIZED 
TO FULLY IONIZED GAS 
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RESUMEN 

Presentamos un atlas 4 de espectros tridimensionales (posicion-posicion- 
velocidad) de la Nebulosa de Orion en tineas de emision opticas de una variedad de 
diferentes etapas de ionization: [O I] 6300 A, [S II] 6716,6731 A, [Nil] 6584 A, [S III] 
6312 A, Ha 6563 A, y [O III] 5007 A. Estas transiciones nos dan information punto a 
punto sobre la estructura fisica y cinematica de la nebulosa a una resolution efectiva 
de 3" x 2" x 10 km s~ , mostrando claramente el comportamiento a gran escala del 
gas ionizado y la presencia de fenomenos localizados tales como flujos colimados 
relacionados a objetos Herbig-Haro. Como un ejemplo de la aplicacion del atlas 
presentamos un analisis estadi'stico de los anchos de las lmeas de Ha, [O III] y [N II], 
que permiten una determination de la temperatura electrdnica media en la nebulosa 
de (9200 ± 400) K. Tambien, a diferencia a trabajos anteriores, encontramos que no 
hay diferencia entre el ensanchamiento no termico de las lmeas de recombination y 
el de las lmeas colisionalcs. 

ABSTRACT 

We present an atlas 4 of three-dimensional (position-position-velocity) spectra 
of the Orion Nebula in optical emission lines from a variety of different ionization 
stages: [O I] 6300 A, [S II] 6716,6731 A, [Nil] 6584 A, [S III] 6312 A, Ha 6563 A, 
and [O III] 5007 A. These transitions provide point to point information about 
the physical structure and kinematics of the nebula at an effective resolution of 
3" x 2" x 10 km s _1 , clearly showing the large scale behavior of the ionized gas 
and the presence of localized phenomena such as Herbig-Haro outflows. As an 
example application of the atlas, we present a statistical analysis of the widths of 
the Ha, [O III], and [N II] lines that permits a determination of the mean electron 
temperature in the nebula of (9200 ± 400) K. We also find, in contradiction to 
previous claims, that the non-thermal line broadening is not significantly different 
between recombination lines and collisional lines. 

Key Words: H II regions — ISM: Herbig-Haro objects — ISM: Herbig-Haro ob- 
jects — ISM: individual: Orion nebula — techniques: spectroscopy — 
turbulence 



1. INTRODUCTION 

The Orion nebula, M42, is the nearest and best 
studied high-mass star-forming region in our galaxy 
(see O'Dcll 2001 and references therein). It is being 
ionized by a compact group of high-mass stars known 
as the Trapezium, which was created in the first 
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of three recent epsiodes of star formation that have 
occured within a region of size one-tenth of a parsec at 
the center of the Orion Nebula Cluster. The other two 
star formation episodes, associated with the BN/KL 
and Orion South regions, are still ongoing and are 
ocurring in dense molecular gas that lies behind the 
ionized nebula. The 07 V star 9 1 Ori C (Simon-Diaz 



et al. 20061 is the most luminous and hottest star 



of the Trapezium and it is responsible for the main 
ionization front. Many morphological features in 
Orion that are apparent in emission lines, such as arcs 
and filaments are associated with collimated gas flows 
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from star-forming activity, such as Herbig-Haro (HH) 
objects (e.g., Reipurth & Bally 2001). Jets, bright 
bars, proplyds and a myriad of other features are 
conspicuous in optical emission lines making of Orion 
an ideal laboratory for studying photoionization and 
hydrodynamic phenomena. Previous high spectral 
resolution studies over restricted regions in Orion 
(e.g.,lCastaneda|1988| [Wen fc Q'Dell||1993| |Esteban 



et al.|1998||Baldwin et al.|2000| ) indicate the presence 



of complex line profiles; the changing complexity of 
the emission line profiles over distinct regions of Orion 
can now be fully appreciated and clearly identified 
with various phenomena over the entire region covered 
by our observations. 

In this paper we present a comprehensive atlas of 
high-resolution echelle slit spectra in multiple emis- 
sion lines, covering a 3' x 5' region in the center of 
the Orion nebula. The atlas combines results from 
two separate datasets, one obtained at Kitt Peak 
National Observatory (KPNO), and the other from 
the Observatorio Astronomico Nacional at San Pedro 
Martir, B.C., Mexico (SPM). Between them, these 
datasets, which have comparable spatial and spectral 
resolution, span a wide range of ionization stages, 
going from [O I] 6300 A, which traces partially ion- 
ized gas at the ionization front, up to [O III] 5007 A, 
which traces the most highly ionized gas in the inte- 
rior of the nebula. Figures and data from the atlas 



are publicly available from the web page http : / /www 
astrosmo . unam . mx/~w . henney/ orionatlas 

Various aspects of the data presented here have 



been discussed in previous papers: Doi et al. (2004 
hereafter Paper I) concentrated on the high-velocity 
features associated with stellar jets; |Henney et al.| 
( 2005 1 analysed the large-scale nebular structure and 



kinematics in terms of numerical simulations of cham- 



pagne flows; Garcfa-Di'az & Henney (2007 hereafter 



Paper I) investigated the global structures found in 



low- velocity emission from low-ionization gas; Hen- 
ney et al. ( |2007 ) analysed the structure and dynamics 
of the low-excitation jet, HH 528. However the full 
set of individual long-slit spectra are shown here for 
the first time and the usefulness and applicability of 
this comprehensive data set goes beyond the analysis 
of the previous papers. As example applications of 
the atlas data, we calculate velocity moment maps 
of the nebula in each emission line and present a 
comprehensive analysis of the line profile statistics. 

The structure of the paper is as follows. In §[2] we 
present the atlas of slit spectra. In §[3] we summarise 
the spatial variations in line profiles across the face 
of the nebula by means of maps of integrated line 
quantities, such as centroid velocity and linewidth. 



TABLE 1 
EMISSION LINES IN THE ATLAS 



Line 



N-S slit positions 
KPNO SPM Total 



[OI] 6300 A 




60 


60 


[S II] 6716,6731 A 


37 


55 


92 


[Nil] 6584 A 


96 




96 


[S III] 6312 A 




60 


60 


Ha 6563 A 


96 




96 


[OIII] 5007 A 


96 




96 



In § [4] we analyse the statistical properties of the 
velocity field in the nebula, including correlations 
between the kinematics in different emission lines. In 
§ [5] we derive the nebular electron temperature from 
the difference in linewidths of Ha and [O III] . In § [6] 
we discuss our results in the light of previous work, 
in §[7] we present our conclusions. 

2. THE ATLAS 

The spectral line data used to construct the atlas 
come from observations at two different observatories, 
KPNO and SPM. The number of north-south ori- 
ented slits from each observatory that contribute to 
the dataset for each emission line is given in Table [l] 

The SPM observations were obtained with the 
Manchester echelle spectrometer (MES; 



Mcaburnl 



et al. 20031 combined with the f/7.8 focus on the 
2.1 m San Pedro Martir telescope of the Universi- 
dad Nacional Autonoma de Mexico, further details 
of the observations are described in Paper 1 1. The 
KPNO observations were obtained with the echelle 
spectrometer attached to the f/8 Cassegrain focus 
of the Mayall 4 m telescope at Kitt Peak National 



Observatory, further details are given in O'Dell et al 



(2001) and Doi et al. (2004) 



The SPM-MES spectra were obtained with a slit 
width of 150 /mi (2.0") for a resolution of 12 km s , 
although a few positions were observed with a nar- 
rower, 70 /im (0.9") slit and a 6 km s _1 , resolution. 
The slit length at SPM is 312". For KPNO the slit 
width is 130 /im (0.8"), the slit length 300" and 
8 km s" 1 resolution. Further details of the data re- 
duction steps are given in | DOH04| and |Paper I| Since 
the original analysis of DOH04 was oriented towards 
fast-moving jets and shocks, the tolerances employed 
for their velocity calibration (~ 4 km s _1 ) did not 
need to be as tight as in |DOH04[ Furthermore, al- 
though |Paper I| employed an E- W slit in order to tie 
together the velocity calibrations of the multiple N-S 
slits, this was not done by |DOH04 Instead, it was 
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assumed that the mean velocity, integrated over each 
slit, was constant as a function of right ascension. 
We have therefore internally re-calibrated each 



DOH04 velocity cube in two ways. First, we have ap- 



plied small rectification corrections in order to remove 
obvious discontinuities between adjacent slits, which 
become apparent when examining the kinematics of 
slow- moving gas. For Ha and [Nil], the affected slits 
are all those lying between x — —55 and x = 47, 
whereas for [O III] no corrections were necessary. A 
common velocity correction, which is linear in dis- 
tance along the slit, is applied to all the affected slits 
in each emission line: SV = a + by (SV is in km s — , 
while x and y are in arcsec with respect to 9 1 Ori C), 
with a = -0.7, b = -0.013 for Ha and a = -0.4, 
b = —0.015 for [Nil]. The velocity corrections to 
individual pixels vary between SV ~ —3 km s" 1 in 
the far south to SV ~ +1 km s _1 in the far north. 
Although this procedure greatly reduced the discon- 
tinuities, it did not remove them completely in the 
case of Ha. 

For the second recalibration, in order to correct 
large-scale trends in the mean velocity as a function 
of right ascension, we have obtained new observations 
in February 2007 through an E-W oriented slit, po- 
sitioned 22"S of 9 l Ori C. These observations were 
obtained with SPM-MES, using a 70 /im slit and 
an instrumental configuration identical to that of 
|Paper I[ except for the filters employed. The mean 
velocities of [Nil], Ha, and [O III] were calculated 
as a function of RA along the slit, and these were 
compared with the mean velocities extracted from 
the |DOH04| velocity cubes for the same slit. For 
[Nil], we found an rms difference < 1 km s _1 , so no 
correction was applied. However, for Ha and [OIII], 
systematic discrepancies of ~ 5 km s _1 are seen in 
the west of the nebula. We have therefore fitted a 
third order polynomial to the velocity differences as 
a function of RA and used it to correct each slit of 
the |DOH04| cube. 

In order to check the calibration of the absolute 
velocity scale for each emission line, we have com- 
pared our results with those from previous studies 
dBalick et al.fl"980] |Castaneda|[l988] |Wen fc O'Dell 
|1993| |Henney fc O 'Dell||1999| [Baldwin et al.pOOOT 
In all cases, we found agreement to within 2 km s~ . 

In summary, we estimate that over most of the 
[S II] , [Nil], Ha, and [O III] maps, the relative veloci- 
ties of each emission line are accurate to ~ 0.5 km s _1 
over relatively short distances (< 50") in the declina- 
tion direction, but only accurate to ~ 1 km s _1 when 
comparing points lying in different slits or widely sep- 
arated in declination. For the [S III] and [O I] lines, 



the low signal-to-noise ratio over much of the map 
means that the accuracy is much worse than this 
for individual spatial pixels, and therefore the same 
velocity precision can only be retrieved by integrating 
the line profile over an extended region and thus los- 
ing spatial resolution. For the absolute value of the 
mean velocity of the different lines, we also estimate 
an uncertainty of 1-2 km s _1 . 

Figures [T}|5] show the individual long-slit spectra 
in the form of calibrated position-velocity (P-V) ar- 
rays. The HST image of Orion at the left of the 
figures shows with a highlighted strip the region cov- 
ered by the sixteen slit positions shown to the right. 
The strip advances to the right (west) in the follow- 
ing rows indicating the corresponding regions for the 
spectra. Figures 1 & 2 correspond to the [S II] 6731 A 
emission line. We do not show here the data for 
the [S II] 6716 A component, although the [S II] line 
ratio has been used by |PaperT] to derive the spatially 
resolved electron density for the whole region. Strips 
A-F in Figures [lland [2]cover the full survey area 
for [SlI]. Figures [3] and [4j (Strips A-C) correspond 
to [OI] 6300 A and finally Figures [Z] (Strip D) and [5] 
(Strips A-B) correspond to [S II] 6731 A. 

The spectra are presented within a heliocentric 
velocity range of —40 to +60 km s _1 in all cases. We 
choose this range because it allows one to appreciate 
the complex structure of the nebula at these relatively 
low velocities, although some high velocity features, 
as those from HH jets fall outside the range displayed, 
but still can be clearly identified. The declination 
scale and right ascension scale set the zero referred 
with respect to 9 1 Ori C. The stellar continua have 
been subtracted from the spectra. Note that the 
mean velocities of the background molecular cloud 
and of the stars in the Orion Nebula Cluster are 
both V© ~ 27 ± 2 km s _1 and that the conversion 
to "local standard of rest" (LSR) velocities is Vlsr = 
V Q - 18.1 km s" 1 . 

3. MAPS OF LINE PROFILE PARAMETERS 

In order to provide a global overview of the nebu- 
lar structure and kinematics, we show in Figures 12 
to [16] maps of derived parameters of the line profile of 
each emission line at each point in the nebula. These 
offer a complementary approach to the isovelocity 



channel maps that have already been published (Doi 



et al |2004||Garcfa-Diaz fc Henney|2007l |Henney et al. 
2007| . The parameters shown are total line surface 
brightness, S, mean velocity, (V), root-mean-square 
(RMS) linewidth, a, peak velocity, V p , and full width 
at half maximum (FWHM) linewidth, W. The first 
three of these quantities can be defined in terms of the 
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Fig. 14. As Fig. [12] but for the peak velocity, V p . 
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first three velocity moments, Mfc, of the line profiles 
I(V): 



V-2 



(1) 



M k = / V«I(V)dV, 

JVi 

such that S = M , (V) = Mx/Mq, and a 2 = 
(M2/M0) — (V) 2 . The limits of integration are cho- 
sen to be Vi = -40 km s _1 and V 2 = +70 km s _1 
for all lines, which represents a compromise between 
being a wide enough range to include most emission 
of interest, but not so wide as to introduce too much 
noise into the maps. This is is a particular issue for a, 
which is very sensitive to noise in the far wings of the 
/(V) profile. The remaining quantities, V p , and W 
are directly measured from the line profile. In order 
to determine Vp to a precision greater than the size 
of the velocity pixels, we use parabolic interpolation 
between the peak pixel and its two neighbours. Simi- 
larly, W is calculated as the difference between the 
velocities of two half-power points, which are found 
from linear interpolation between adjacent pixels. In 
the case that the profile has more than two half-power 
points, W is defined as the difference between the 
highest and lowest of these. 

Figure [12] shows the reconstructed surface bright- 
ness maps in the six lines. Comparison with pub- 
lished direct imaging (e.g., Pogge et al.|1992 l shows 
generally good agreement, although some artefacts 
remain from the process of combining the individual 
spectra. These can easily be identified as brightness 
discontinuities across vertical lines in the maps. The 
most serious of these artefacts are in the [S II] map 
at x = +10, x — +45, in the [Nil] and Ha maps 
at x = +45, in the [S III] map at x — 0, and in the 
[O III] map at x — — 5. 

The mean and peak velocities are shown in Fig- 
ures |T3l and [141 where it can be seen that the two 



quantities are generally well-correlated for a given 
line, although there are some important differences. 
The mean velocity varies smoothly across the face of 
the nebula, whereas the peak velocity shows frequent 
discontinuities, which are caused by one particular 
line component becoming brighter than another. The 
mean velocity is also more sensitive to high velocity 
components, such as HH objects. This is particularly 
notable in the case of [O III], where the HH objects 
do not show up at all in the V p image. Since [S III] 
and [O I] are very much weaker than the other lines, 
the velocity maps are rather noisy, especially in the 
faint regions of the nebula. We have therefore masked 
out from the velocity maps all regions with a surface 
brightness below a certain threshold, which appear 
as solid gray areas in the figures. 

The non-thermal contributions to the RMS and 



FWHM line widths are shown in Figure [15] and 16 
These are calculated by subtracting in quadrature 
the instrumental and thermal widths from the raw 
measured widths, assuming a gas temperature of 
9100 K (see §[5|. For a Gaussian profile, the rela- 
tionship between the FWHM and RMS widths is 
W = (81n2) 1 / 2 cr ~ 2.355cr. We therefore multiply 
the values of a by 2.355 to allow direct comparison 
between the a and W maps. There is generally a 
much greater difference between the a and W maps 
of a given line than between the maps of (V) and 
V p . The a maps are very much dominated by high- 
velocity flows, which are much less apparent in the 
W maps. As with the V p maps, the W maps show 
discontinuities between patches of the nebula with 
high and low line widths. An exception is the Ha 
line, where the increased thermal broadening leads 
to a more Gaussian-like line profile and a closer re- 
semblance between the W and a maps. As with the 
velocity maps, we have masked out regions of low 
signal-to-noise ratio for the [O I] and [S III] lines. 

An empirical description of the nebular structure 
and kinematics was given in § 3 of Paper I and will not 
be repeated here. The reader is referred in particular 
to Fig. 7 of Paper I for orientation in identifying 
nebular features mentioned in the remainder of the 
paper. 

4. ONE-POINT VELOCITY STATISTICS 
In Figures [T7| and fl8] we show joint distributions 
of the mean and peak velocities, respectively, for se- 
lected pairs of emission lines, which illustrate the 
correlations that exist between the kinematics of the 
different emission zones of the nebula. Figure \l9\ 
shows the joint distributions of the raw velocity vari- 
ance, a 2 , which is the square of the RMS velocity 



dispersion, and Figure 20 shows distributions of the 
non-thermal FWHM linewidth, W . The values of 
W have been approximately corrected for thermal, 
instrumental and fine-structure broadening in the 
same way as in Figure 16 whereas the values of a 2 
in Figure |19| are uncorrected and measured directly 
from the data. 

The statistics for individual lines are summarised 
in the upper part of Table [2] which shows the mean 
and standard deviation (s.d.) of each quantity. These 
are calculated with an equal weighting of each pixel 
in the map after masking out any points lying out- 



side the ranges shown in Figures 17 20 For [S III] 
and [O I] , any pixel with surface brightness below the 

TBI were also masked 



thresholds shown in Figures [13} - 
out. The lower part of Table [2 show the mean and 
s.d. of the differences between the kinematic quanti- 
ties for the line pairs shown in the figures, together 
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Fig. 17. Joint distribution of mean velocities, (V), of different emission lines, with inset box showing a histogram of 
differences in (V). In each case, the straight line shows the case of equal velocities in the two lines, (a) {V)ua versus 
(V)[onrj- (b) (V)[ N ii] versus (F)[oiii] (c) (V)[sirj versus (V)[ S ni]. (d) (V) [0 i] versus (V)[ S ii] 
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Fig. 18. Joint distribution of peak velocities, V p , of different emission lines, with inset box showing a histogram of 
differences in V p . In each case, the straight line shows the case of equal velocities in the two lines, (a) Vp(Ha) versus 
VpQOIII]). (6) Vp([NII]) versus Vp([OIII]) (c) Vp([SII]) versus VpQSIII]). (d) V p ([Ol]) versus V P ([S II]) 
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Fig. 19. Joint distribution of raw measured mean square velocity widths, a 2 , of different emission lines, with inset box 
showing a histogram of the differences in a 2 , (a) versus cr 2 OIU y (b) <7p NII ] versus cr 2 olll y (c) crp SII ] versus of SIII ]. (d) 
°fo i] versus a 2 s n j . The diagonal straight line in each panel shows the case of equal variances in the two lines, except for 
in panel a, where it shows the relation 4„ = ofoin] + 82.2. 
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Fig. 20. Joint distribution of corrected full-widths at half maximum, W, of different emission lines, with inset box 
showing a histogram of differences in W. The thermal and instrumental widths have been subtracted in quadrature from 
the raw measured values. In each case, the straight line shows the case of equal widths in the two lines, (a) Whc versus 
W[om]- (b) W [NU] versus W [0 ni] (c) W [S n] versus W [S ni]. (d) W [Q i] versus W [su] 
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TABLE 2 

SUMMARY OF VELOCITY STATISTICS 





line 


(V) 




a 2 


W 




[OIII] 


16.3 ±2.8 


15.5 ±4.0 


113 ±29 


14.7 ±5.1 




Ha 


16.8 ±3.0 


16.2 ±4.1 


195 ± 28 


18.6 ±5.1 




[Sin] 


19.7 ±3.0 


19.2 ±4.0 


99 ±46 


14.6 ±4.4 




[Nil] 


20.5 ±2.9 


21.1 ±3.7 


106 ± 24 


15.4 ±3.7 




rc ttI 

[b iij 


01 O I O 4 

zl.z ± z.4 


21.0 ± 3.2 


11/) I 1 

114 ± zl 


ICO _L O C 

15.2 ± 3.5 




IO Tl 
IV fl 






1UU —1— U<J 


11.5 ± 2.6 


Ha — 


[0 III] 


0.5 ± 1.4 


0.7 ± 3.2 


82.2 ± 11.8 


3.9 ± 3.0 




r 


0.889 


0.700 


0.917 


0.825 


[N II] — 


[OIII] 


4.2 ±2.8 


5.6 ±4.5 


-7.0 ±36.8 


0.7 ±6.5 




r 


0.511 


0.274 


0.287 


-0.067 


[S II] - 


[S III] 


1.6 ±2.5 


2.6 ±3.3 


18.0 ±46.8 


0.2 ±5.6 




r 


0.617 


0.588 


0.229 


0.065 


[01]- 


- [S II] 


4.4 ±3.0 


4.5 ±2.7 


-7.6 ±64.8 


-3.7 ±4.2 




r 


0.507 


0.580 


0.175 


0.078 



with the linear correlation coefficient, r, of the joint 
distributions. In cases where the quantities are well 
correlated (r > 0.5), the s.d. of the difference can 
be significantly less than the s.d.'s of the individual 
values. 

The general trend of increasing blueshifts with 
increasing degree of ionization is readily apparent 
from the table, which is arranged from highest to 
lowest ionization. The mean and peak velocities 
always show a positive correlation between different 
lines. The correlation is highest (r ~ 0.9) between 
Ha and [O III] , which is not surprising since there is a 
large overlap of the emission zones of these two lines. 
However, even pairs of lines that should not greatly 
overlap (such as [N II] - [O III] and [S II] - [S III]) also 
show significant correlation (r > 0.5). 

The s.d. of 1.4 km s _1 for (V)h q - (^)[om] is a 
hard upper limit on the uncorrelated average wave- 
length calibration errors within the individual lines. 
Given that (U)[nh] — (V)[om] has a much larger 
s.d. (2.8 km s _1 ) and that some fraction of the Ha 
emission comes from the [N II] zone, it follows that at 
least some of the scatter in (V)n a — (Varolii] must 
be real. This is consistent with our estimate above in 
§ [2] that the relative velocity calibration within each 
emission line is good to ~ 1 km s _1 . Our worst-case 
estimate for the uncertainty in the absolute velocity 
calibration for each line is 2 km s _1 and this trans- 
lates into uncertainties in the average values of the 
differences in mean and peak velocities. This issue 



is returned to in § [5] below, where we argue that the 
absolute velocity calibration must be better than this 
conservative estimate. 

An extremely tight correlation (r ~ 0.9) is seen 
between the mean square velocity widths of Ha and 
[O III] , whereas no other pair of lines show a signif- 
icant correlation in this quantity. The large offset 
(0H a — ofoiiil) = (82.2 ±18.2) km 2 s~ 2 is due mainly 
to the extra thermal broadening of Ha, as is discussed 
further in §[5] below. Other pairs of lines show much 
smaller offsets in cr 2 , with no obvious trends with 
ionization. The s.d. of a 1 is very large in the case of 
[O I] and [S III] because of the low signal-noise ratio 
of these weak lines. 

A similar picture is seen for the corrected full 
width half maxima, W: only in the case of Ha — 
[O III] does a correlation exist. Note that even though 
the effects of thermal and fine structure broadening 
have been subtracted in quadrature to obtain W, a 
residual difference remains between Ha and [O III] 
of ~ 4 km s _1 . We believe that this is because the 
FWHM is such a mathematically ill-behaved quantity 
that quadrature subtraction cannot be expected to 
work reliably in the presence of strongly non-Gaussian 
lineshapes. For example, consider a line that consists 
of two components, A and B, separated by 5v and 
each of width w, where the intensity of B is less 
than half that of A. The measured FWHM will not 
be affected at all by the presence of component B 
unless the components are blended: Sv < w. The 
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larger thermal broadening of Ha gives a larger w, so 
that there will be cases where the FWHM of Ha can 
"see" component B, but the FWHM of [O III] cannot. 
Evidence for this can be seen in Figure |16| where 
the [O III] map shows multiple dark patches of low 
FWHM, which are absent from the Ha map. 

All the metal lines show very similar mean values 
of (W) ~ 15 km s _1 , except for [01], which is sig- 
nificantly lower, (W) ~ 12 km s _1 . The FWHM is 
affected by noise to a much smaller extent than a 1 
and the metal lines all have low thermal widths, so 
should not be affected by the problem discussed in 
the previous paragraph. 

5. DERIVATION OF THE MEAN ELECTRON 
TEMPERATURE FROM LINEWIDTHS 

In this section, we use the statistics of the ob- 
served line profiles of [Nil], Ha, and [O III] (§(4j| to 
derive the mean electron temperature in the nebula 
and to show that the "non-thermal" component to 
the line broadening is not significantly larger in Ha 
than in the metal lines. The method we employ is 



essentially an old one (Courtes et al. 1968 Dyson 



fc Meaburn||1971 1 , which estimates the gas temper 



ature by comparing the width of a metal line with 
that of a hydrogen line, under the assumption that 
the excess width of the hydrogen line is principally 
due to its greater thermal broadening on account 
of its lower atomic weight. A problem with many 
applications of this method has been that the metal 
line and hydrogen line (usually [N II] and Ha) are 
not necessarily emitted by the same volume of gas. 
Indeed, previous applications of this method to the 
Orion nebula ( Dopita et al.|1973 Gibbons|1976 1 have 
shown that the H/3-[0 III] pair gives results that are 
very different from the Ha- [Nil] pair and that are 
seemingly more reliable. Our own version of the 
method improves on previous versions by using all 
three emission lines [Nil], Ha, and [O III] in order 
to control for the imperfect overlap of the emission 
regions. 

5.1. Components of the observed line widths 

For each emission line, the observed mean square 
velocity width at a certain point on the nebula can 
be broken down schematically into four components: 



'th 



9 



(2) 



The first term in this equation is the thermal Doppler 
width, cr t 2 h = 82.5(T 4 /A) km 2 s~ 2 , where T 4 = 
T/10 4 K and A is the atomic weight of the atom 
or ion. The second term is the fine structure broad- 
ening, which is important for hydrogen and helium 



recombination lines, but not for metal lines such as 
[O III] and [N II] . In Appendix |Bj we show that its 
value for Ha is af s = 10.233 km 2 s 2 . The third term 
is the instrumental width, which for the KPNO ob- 
servations is approximately of ns = 11.5 km 2 s , and 
which is the same for the three lines. The fourth term, 
cr 2 t , is the "non-thermal" contribution, which is a 
catch-all term for any other broadening process. This 
might include Doppler broadening due to velocity 
gradients or turbulence within the emitting regions, 
broadening due to scattering by dust particles, or any 
other process. 

Applying the above equation to differences in a 2 
measured in § H] gives 



2 2 
a Ha - a [0 III] 



77.34T 4 + 10.233 



and 



°fNii] _ <J ronil - 0.74T 4 



[o III] 



CT nt,Ha _ CT nt,[OIII] (3) 



' a nt,[NII] _ CT nt,[OIII]- (4) 



Note that the instrumental contribution cancels from 
these equations. The simplest way of extracting the 
temperature, and what has been done in the past, 
is then to assume that the difference in non-thermal 
widths is negligible, so that T 4 can be calculated 
directly from the measured o\ a — er 2 Ilr i as 



rpl 



'Ha 



a [OIII] 



10.233) /77.34 (5) 



However, by including the data on the [N II] line, it 
is possible to derive a lower bound on the value of 
and so arrive at a more reliable 



'nt.Hc 



"ntJOIII] 



temperature determination than the naive estimate 

5.2. Correction for differences in temperature and 
non-thermal widths between high and low 
ionization zones 

We adopt a simplified two-zone model for the 
emission structure along a given line of sight through 
the nebula. The first zone, A, contributes all the 
[O III] emission plus a fraction, /, of the Ha emission. 
The second zone, B, contributes all the [N II] emission 
plus the remaining fraction, 1 — / , of the Ha emission. 
We further assume that the non-thermal broadening 
is the same for all emission lines that arise from the 
same volume (this assumption should hold true for 
turbulent or other kinematic broadening, an issue 
we will return to later). A fundamental property of 
the moments of the line profile given in equation ([lj 
is that they are all additive, with consequences that 
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are elucidated in Appendix [A] From equation ( A9 1 
it follows that 



2 

a nt.Ha 



r nt,[OIII] — (1 _ /)( CT nt,[NII] ~ °"nt,[OIII]) 

+ /(1-/)(W[nii]-W[oiii]) 2 - (6) 



Note that the above equation is only valid if the 
effects of dust extinction are the same for zones A 
and B, but this is probably a good approximation 
since 10'Dell & Yusef-Zadehl (12000b showed that the 



majority of the dust extinction in Orion arises in the 
neutral veil, in the foreground of the ionized gas. 

We also allow the two emission zones to have 
different temperatures, T A and Tb, such that the 
mean Ha temperature is T = fT A + (1 — f)T&. We 
can now combine equation ([6]) with equations ^ 
and Q, modified to account for the temperature 
difference between the zones, to obtain 



T 4 = rjTi + e - C - & 



(7) 



where 



ry= (1 - 0.00957 (l-/))~\ (8) 
e = 0.00957 77/(1 - f)(T B - T A )/10 4 K, (9) 
C = 0.01293 7,(1 - /)(af NII] - af om] ), (10) 

£ = 0.01293 77/(1 - /)«V)[nii] - (V) [0 m]) 2 - (11) 

5.3. Techniques for estimating the ionization 
factor f 

Equation Q gives the temperature at any point in 
the nebula in terms of observable quantities and the 
factor /, which can be estimated in two ways: from 
surface brightness ratios, or from ratios of differences 
in mean velocity. We shall describe each of these 
techniques in turn. 

We define for each emission line a per-zonc 
mean effective emission coefficient, /3, with units 
erg cm 3 s _1 sr _1 , such that the intrinsic surface 
brightness of a zone is 5* = (3E, where E — J n p n e dz 
is the emission measure of the zone. Since [N II] 
6584 A and Ha 6563 A are so close in wavelength, 
they are both affected in the same way by dust ex- 
tinction, so the observed ratio of their surface bright- 
nesses is equal to the intrinsic value. Therefore, given 
that [Nil] by construction arises only in Zone B, we 
have 



R 



[Nil] 



S([NII]) = /3b([NII])£ b 
S(Ha) (3 A (Ha)E A + (3 B (Ra)E B 
^bQNII]) 



ii-fV 



/3b (Ha) 



(12) 



where the last equality follows from the definition 
of /. Hence, if the emission coefficients /3b ([Nil]) 
and Pb (Ha) are known, then / can be calculated 
from the observed surface brightness ratio as / = 
1 - i?[Nii]/3B(Ha)//3 B ([NH]). To make use of this 
equation, it is first necessary to photometrically cali- 
brate the surface brightness maps and we have done 
this by comparison with published spectrophotomc- 



try of the region 


Pogge et al. 


1992 


Baldwin et al. 


2000 


O'Dell et al. 


20031. It is then necessary to esti- 



mate the ratio of emission coefficients for conditions 
typical of Zone B: T ~ 10 4 K, n c ~ 5000, which we 



have done using the Cloudy plasma code ( Ferland 



20001, assuming an N/H abundance ratio of 6 x 10 



finding that the ratio is roughly unity. This is con- 



sistent with the N + /N abundance found by Esteban 



et al. ( 1998 1 



Assuming that /3 B (Ha)//3 B ([N II]) = 1, we find 
a mean value of (/) = 0.76 for the entire area cov- 
ered by our observations. The upper-left panel of 
Figure [21] shows the map of / calculated by this 
technique, which varies between 0.26 and 0.91. In 
reality, (3b (Ha) //3b ([N II]) will vary with position 
due to variations in electron density ([Nil] is af- 
fected by collisional deexcitation but Ha is not). 
A better way of calculating / would be to use the 
[OlII]/H/3 surface brightness ratio, since the criti- 
cal density for [O III] is higher than that typically 
found in the nebula (similar arguments to those 
above give / = R [OIII] f3 A (H.f3)//3 A ([0 III]), where 
R [om] = 5([OlII])/5(H/3)). However, our dataset 
does not include H/3 and the [O III] /Ha ratio would 
have to be corrected for the highly variable fore- 
ground extinction ( |Q'Dell fc Yusef-Zadeh||2000[ ) . A 
comparison of the [N II] /Ha and [O III] /H/3 ratios 
over a limited area (using data from Figs. 11 and 12 
of |Q 'Dell et~aT|2003| gives an rms scatter of ±0.05 
between the values of / derived from the n] and 
-R[Oiii] techniques, which is roughly half of the to- 
tal rms variation in / from our maps. The -R[om]~ 
derived values of / have a greater range than the 
i?[ N u]-derived values. This is to be expected since, 
all else being equal, higher density regions will have 
lower ionization. 

A completely independent way of estimating /, 
and one that does not depend on knowing the emis- 
sion coefficients, is to use the observed mean velocities 
of the three lines. A corollary of the additive property 
of the line profile velocity moments (Appendix |A| 
is that in the two-zone model we have from equa- 



tion (A8l that 



(V)b q = f(V) l0 ui] + (1 - f)(V) [NU ], (13) 
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0++ fraction, / Correction factor Electron temperature, 'l' c /10 4 K 



0.25 0,5 0.75 1 -0.25 0.25 0.5 0.75 1 1.25 
^M^t^l * I I ^^^^M I I II 




-100 -50 50 -100 -50 50 -100 -50 50 
x offset from 8 l C (") x offset from & 1 C (") x offset from 6 i C (") 



Fig. 21. Steps in the construction of the temperature map. Left column: Two independent estimates of the fraction, /, 
of the Ha emission that comes from the ++ zone: from the [Nil] /Ha surface brightness ratio (upper panel), and from 
the differences in mean velocities of [Nil], Ha, [O III] (lower panel). Center column: Correction factors to the "naive" 
temperature map due to differences in the kinematics between the N + and ++ zones: mean velocity differences, £ 
(upper panel) and velocity dispersion differences, £ (lower panel). Right column: "naive" temperature estimate (upper 
panel) from the difference in velocity dispersion of [OIII] and Ha and final corrected temperature estimate (lower panel). 
The horizontal bands and vertical streaks seen in some panels are instrumental artifacts. 
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so that / can be estimated as 



/ = 1 



+ 




2 4 6 

<V s >[Nn] - (V B ){am], 



km f 



Fig. 22. Joint distribution of differences in mean velocities: 
{V)u a - (F}[oiii] versus {V)[ N ii] - (V)[oin]. The values 
of (V}Ha — {V)[o in] have been shifted by 0.5 km s _1 so as 
to make the distribution symmetric about a line through 
the origin (see discusion in text). The dashed lines show 
the boundaries of the two triangular regions (upper right 
and lower left) that give physically possible values of / in 
the two-zone model. The solid line shows the relationship 
(V) Ha - (V) l0 m] = (1- (/»((V) [NII] - (V>[oni]), where 
(/) is determined from the -R[nii] method. The solid gray 
circle indicates the effect of ±1 km s _1 uncertainties in 
the velocities. 



<V)Ha-<V)[QIII] 
00 [Nil] - (^)[OIII] 



(14) 



Despite its advantages over the surface brightness 
method, this technique has two serious limitations. 
First, equation ( 14 ) has a singularity when (V) r N n] = 
(V) [o in] • If the two-zone model were strictly valid, 
then any pixel with (U)[nii] = (U)[oni] would also 
have (U)hq = (^0[oiii]j so the singularity would be 
benign. However, in reality the two-zone model is 
only an approximation and, additionally, the observed 
values of (V) will always include a certain measure- 
ment error. As a result, applying equation (14 1 to 



real data will inevitably result in discontinuities in 
the / map and unphysical values (/ < or / > 1) 
for some pixels. Second, the observed values of 
(V)hu — (V)[oui] are rather small (typically —2 to 
+3 km s _1 , see Fig. 17 1), and are of the same or- 



der as the systematic uncertainties in the wavelength 
calibration. 

Given these caveats, it is perhaps surprising that 
this method works at all. However, the average slope 
of the relationship between (V)n a — (^0[om] an( f 
(^)[nii] ~~ (^0[om] should be independent of any 
inaccuracies in the absolute velocity scales of the 



individual lines, and by equation ( 14 1 should be equal 
to 1 — (/). This is illustrated in Figure 22 where 
it can be seen that the value (/) = 0.76, calculated 
from the jji method above and shown by the solid 
line, is indeed close to the mean trend in the velocity 
differences. Furthermore, to the extent that the two- 
zone approximation is valid, the symmetry axis of the 
joint distribution of velocity differences ought to pass 
through the origin. This is only true when we shift 
all our (V)n a — (V)[om] values by +0.5 km s 1 and 
this is what is done in the figure. Such a small shift is 
well within the estimated uncertainties in the velocity 
zero point of each line. Of course, there may also 
be a similarly small constant error in the (V)[Nn] — 
(V}[oiii] values, but we have no information that 
would allow us to determine the correction in this 
case. 

In order to correspond to a value of / between 
and 1, the velocity differences must lie in one of 
the two triangular areas of the graph in Figure [22] 
bounded by dashed lines. It can be seen from the 
graph that the majority of pixels do indeed lie in 
these areas, mostly on the positive side of the origin. 
The small fraction of pixels lying in unphysical areas 
of the graph can be explained by residual errors in 
the velocity determination of order 1 km s _1 (indi- 
cated by the solid gray disk). The lower left panel 
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TABLE 3 

SUMMARY OF TEMPERATURE 
DETERMINATIONS 







(2*) a 


t 2 


AW 


naive 


0.930 


0.942 


0.025 


178432 


corrected 


0.913 


0.916 


0.014 


178432 


low a 


0.916 


0.919 


0.012 


173180 


high a 


0.824 


0.792 


0.104 


5252 



of Figure [21] shows the map of / calculated using 
this velocity difference technique. For much of the 
map, the agreement with the surface brightness ra- 
tio technique (upper left panel) is rather good, and 
the low ionization regions associated with the linear 
bright bars are readily apparent in both maps. The 
unphysical values of / show up as solid black or white 
regions in the map. The fact that these regions occur 
predominantly at the edges of the map, with white 
regions (/ > 1) concentrated in the West and black 
regions (/ < 0) in the East, suggests that the resid- 
ual velocity errors are largely systematic, rather than 
random. 

5.4. Corrected value of the electron temperature 

We are finally in a position to calculate the 
correction factors, e, f, and £ needed to estimate 
the electron temperature. The mean values are 
(e) = 2.9 x 10~ 4 , (C) = -0.0218 and (£) = 0.0417, 
so that the corrected mean temperature from equa- 
tion Q is (J4) = 0.913, whereas the naive estimate 
would have been (T4) = 0.930 (note that all these 
averages are calculated with equal weighting for each 
pixel). Maps of the spatial variation of the correction 
factors £ and £ are shown in the central panels of 
Figure [2l] while the right-hand panels show maps of 
the naive and corrected temperatures. The correction 
factor e is uniformly very small and will not be con- 
sidered further. The effect of the corrections on the 
mean temperature is rather small, merely lowering it 
by 200 K, but, as can be seen from the maps, the cor- 
rection factors can be much larger in localized regions 
(> 2000 K). Some of the variation seen in the uncor- 
rected T4 map is removed by the correction factors, 
with the most notable example being at the bright 
bar, which is seen as a dark diagonal band in the T4 
map, but not at all in the T4 map, principally due to 
£. Another example is the "Red Fan" ( Garcia-Dfaz 



fc Henney||2007| , seen as a light feature in the lower 
right of the T' A map, but which is effectively removed 
by the £ correction (the thin vertical streaks seen 
superposed on this feature are observational artifacts, 



which can also be seen in the <Jn a map of Fig. 15 1. 
Other striking features seen in the temperature maps 
are those associated with the high-velocity jets em- 



anating from the Orion South region (Henney et al 



2007[ ) : HH 202 (northeast), HH 203/204 (southwest), 
and HH 529/269 (east- west). These features also 
tend to have high absolute values of £ or £, but the 
application of the correction factors is not successful 
in removing them from the temperature map. De- 
spite this, we do not believe that they correspond 
to real temperature variations. Since all these fea- 
tures have large rms velocity widths (Fig. 



151, we 



have masked out of the corrected map all pixels with 
"foini or a [Nii] g rea ter than 180 km 2 s 2 (visible as 
small solid gray areas in the figure). 

In order to quantify the variation in the naive and 
corrected temperature maps, we calculate the plane- 
of-sky Peimbert t 2 parameter ( Peimbert|l967 O'Dell 



et al 



2003| ): t 2 A = ((T-(T) E )*)J(T) B . The aver 
ages m this expression, (■ ■ -} E , should theoretically 
be weighted by the emission measure of each pixel, 
but we approximate this by weighting by the Ha sur- 
face brightness: (• • • ) s . The results are summarised 
in Table [3] It can be seen that the Ha-weightcd 
averages, (24)3, are hardly any different from the 
uniform- weighted averages, (X4). The fractional vari- 
ance, t\, is significantly reduced after applying the 
temperature corrections (second row of table). A 
further small reduction in t\ occurs (third table row) 
when the high-er pixels are eliminated as discussed 
above. The high-er pixels themselves (fourth table 
row) show a significantly lower mean temperature 
and a very large t\. 

6. DISCUSSION 

6.1. Comparison with previous mean temperature 
determinations 

Previous determinations of the mean value of 
the electron temperature in the Orion Nebula have 



been made using various methods (see Osterbrock & 
Ferland 2006 Chapter 5 for an overview). In this 
section, we compare our own results (§ 5.4 ) with some 
recent examples of these. 



O'Dell et al. (20031 used HST emission line imag- 
ing m the [O III] 4363 A and 5007 A lines to determine 
2~[Oin] = 8890 ± 496 K, where the scatter represents 
real variations, rather than observational uncertain- 
ties. Before comparing with our own result, this 
temperature needs to be corrected in two ways. First, 
temperature fluctuations within the [O III]-emitting 
region will lead to this method overestimating the 
real temperature there, T(0 ++ ). Second, T(0 ++ ) 
(which corresponds to Ta from § 5.2 above) is an un- 
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derestimate of the mean temperature of the ionized 
gas, T(H + ), since it does not include the low ioniza- 
tion zone along each line of sight (Zone B), which 
tends to have a higher temperature. Taking estimates 
of the magnitudes of these corrections from |Esteban| 
et al.|(|1998| and|Rubin et al.| ( |2003[ ), we find that the 
O'Dell et al.| fl2003|) resultslmply T(0++) = 8400 K 



and T(H + ) = 9100 K, which is very similar to our 
result of 9190 K ( Table [|). 

Electron temperatures have also been derived 
from the brightness relative to the continuum of radio 



recombination lines ( 


Wilson & Pauls||1984 


Wilson & 


Jaeger||1987 


Wilson et al.||1997l, consistently finding 



values of around 8500 K. Yet another technique is 
to measure the strength of the Balmer dis c ontinu ity 



at 3646 A, which was used by Liu et al. (19951 to 
measure a mean value of T ~ 9000 K, whereas |Es-| 
teban et all (119981) used the same technique to find 



T ~ 8600 K. In the presence of temperature fluctu- 
ations, these values are underestimates and should 
be increased by ~ 300 K. However, the experimen- 
tal uncertainties of these results are probably higher 
than for the other methods and furthermore they 
only correspond to restricted regions of the nebula. 

In summary, once temperature fluctuations are ac- 
counted for, there is a very good agreement between 
the three independent optical determinations of the 
mean electron temperature (line ratios, Balmer dis- 
continuity, and line widths), giving an average value 
of ~ 9100 ± 100 K. The radio-derived temperatures 
are significantly lower by about 600 K and we have 
no satisfactory explanation for this discrepancy. 

6.2. Plane-of-sky variation in the electron 
temperature 

The point-to-point variations in electron temper- 
ature across the face of the nebula have been studied 



in detail by two recent papers (O'Dell et al. 2003 
Rubin et al.||2003 1 based on optical line ratios from 
HST observations. The values obtained for t\{0 ++ ) 
for differ ent regions of the nebula range from 0.005 

find- 



to 0.016. 



Rubin et al. 



also calculate i^(N ++ ), 



ing broadly similar values. After correcting for the 
effects of noise, |Q'Dell et al.| find a value for their 
entire surveyed region (roughly one third the area of 
our own dataset) of t\{0 ++ ) = 0.0079. Given that 
our own value of t\ = 0.012 will also be somewhat 
affected by noise, the agreement may at first seem 
satisfactory. However, the spatial resolution of our 
own observations is at least an order of magnitude 
lower than these two studies, which find that power 
spectrum of the fluctuations peaks at sub-arcsecond 
scales. Such fluctuations would be entirely smoothed 



out by atmospheric seeing in our own dataset, so one 
would expect our results to show a much lower value 
cfii. 

An analysis of the temperature variations along 
the slits used in the |Rubin et al.| study supports 
this point of view. Their Slit 5 is the only one that 
does not cross high-velocity jet features, and we find 
t\ = 0.002 for this slit, which is very much lower than 
their value of 0.018, as is to be expected from our 
lower spatial resolution. It therefore seems likely that 
much of the variation visible in even our corrected 



temperature map (lower right panel of Fig. 21 1 is 



spurious, and that t\ = 0.002 is a more accurate 
estimate of the plane-of-sky temperature fluctuations 
at scales > 1 arcsec. This is also consistent with 
calculations of the variation in three small 40" x 40" 
patches of our map, which were selected for their high 
signal-noise ratio and avoidance of high- velocity flows. 
The resulting values are t\ = 0.002, 0.004, and 0.002 
for regions centered at (—30, +60), (—30,-100), and 
(—80,-50), respectively. 

6.3. The nature of the non-thermal line broadening 

A fundamental assumption of our method of de- 
terming the temperature is that the non-thermal 
broadening depends only on the emitting volume 
(see § 5.2 1, and does not vary systematically between 
collisionally excited and recombination lines. This 
is necessary so that (7 2 t H can be determined from 
the kinematics of the [Nil] and [O III] lines, as in 
equation (J6|. The fact that we determine a mean 
nebular temperature that closely agrees with that 



from other techniques (§ 6.1 1 is a strong argument 
that this assumption is valid. 

However, this is in conflict with the result of 



O'Dell et al. (20031, who found non-thermal widths 
(FWHM) of ~ 19 km s _1 for recombination lines 
(Ha, H/3, and He I 5876 A), but only ~ 11 km s _1 
for collisional lines ([Ol], [S II], [Nil], [S III] , [OIII]). 
Taken at face value, these widths would imply a large 
difference of (19 2 - ll 2 )/(81n2) = 43 km 2 s~ 2 be- 
tween cr 2 t Ha and <7 2 t r nl j , meaning that over half 



of the difference er 2 ^ 



u [o in; 
In this case 



(Fig. 



non-thermal effects 
temperature would only be 



19 a) is due to 



the mean electron 
4000 K, in stark con- 
trast to all other temperature determinations for the 
nebula. 

Instead, we suggest that the difference in non- 
thermal widths found by |Q'Dell et al.| is an arti- 
fact of their methodology, which consisted in fitting 
Gaussian components to the line shapes. The basic 
problem with this is that the extra thermal broad- 
ening of the H and He lines tends to blend the fine 
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details of the line shape, allowing a satisfactory fit to 
be obtained using fewer Gaussian components than 
are necessary for the metal lines. These fewer com- 
ponents would then naturally tend to have larger 
individual widths (even after correcting for thermal 



broadening), which may explain the O'Dell et al. 
result. Our own method of calculating the RMS 
linewidths does not suffer from this problem, since 
it avoids any subjective decision about how many 
components to fit. Indeed, an application of our 
methodology to the same data as used by |Q'Dell 



let al.l shows no evidence for extra non-thermal broad- 



ening of the recombination lines (Fig. 4 of Henney 
fc Q'Dell|[l999| . A further problem with the |Q 'Dell 



et al. ( 2003 1 study is that no correction was made for 



the H fine-structure broadening (although the larger 
He I fine structure was accounted for). 

A third methodology was employed by |Gibbons] 
who convolved the observed [O III] profile 



(1976) 



with Gaussians of different widths and determined 
which width gave the best match to the observed H/3 
profile. This technique gave temperatures that are 
very similar to those derived here, which is further 
evidence that the IQ'Dell et al.l widths cannot be 
correct. 

One promising class of explanations for the non- 
thermal linewidths is that they are kinematic in na- 
ture, caused by velocity structure in the ionized gas. 
This may be due to velocity gradients in ordered, 
large-scale champagne flows (e.g. 



Yorke et al. 1984 Henney et al. 



Tenorio-Tagle 1979 



2005), or to multi- 



Henney et al. (2005) 



2006 



scale, turbulent motions (e.g. |Mellema et al 
However, 

showed in their § 6.1.2, velocity gradients in an or- 
dered flow are incapable of explaining the observed 
linewidths. Furthermore, the disordered appearance 
of the maps of the line profile parameters (Figs. 13 

In 



to 



16 1 would favor the turbulent explanation. 



the radiation-hydrodynamical simulations of Mellema| 
(2006), the turbulence in the ionized gas is 



et al. 



driven primarily by photoevaporation flows from 
dense concentrations of molecular gas, which are 
themselves the product of supersonic turbulence in 
the molecular cloud from which the high-mass star 
formed. The internal velocity dispersion of the H II 
region remains roughly sonic for the full duration of 
the simulations (~ 0.5 Myr), in agreement with the 
observed linewidths in Orion. In addition, it is quite 
possible that other mechanisms may be contributing 
to the turbulent driving, such as ionization front in- 
stabilities ( Garcia-Segura k. Franco|[l996 Williams 



1999 Whalen & Norman 2007) or supersonic jets 



7. CONCLUSIONS 

We have presented an atlas of spatially resolved 
high resolution emission line spectra of the Orion 
Nebula, covering a wide range of degrees of ioniza- 
tion. This has allowed us to investigate the global 
kinematics of the nebular gas in unprecedented depth. 
We have studied the statistical correlations between 
the line profile parameters of the different emission 
lines. Our principal conclusions are as follows: 

1. The mean velocities of different emission lines 
are always positively correlated. The correlation 
is greatest when there is considerable spatial 
overlap between the emitting volumes of the two 
lines, but correlation also exists between lines 
from non-overlapping volumes. This is indicative 
of the existence of an ordered mean flow that 
connects the spatially separate regions. 

2. There is no significant correlation, on the other 
hand, between the linewidths of emission lines 
from non-overlapping volumes, such as [N II] 
6584 A and [O III] 5007 A. This is consistent 
with kinematic broadening of an essentially tur- 
bulent nature. The kinetic energies associated 
with both ordered and unordered motions are of 
the same order as the thermal energy. 

3. The mean electron temperature in the nebula is 
determined to be 9190 K from comparison of the 
widths of Ha and [O III] lines, after correction 
for differences in non-thermal widths between 
the different emission zones. Although apparent 
spatial variations in temperature are detected 
{t\ = 0.012), most are probably due to uncer- 
tainties in these corrections, rather than real 
temperature variations. In regions of the nebula 
where these correction terms are small, the spa- 
tial variations are much smaller (t\ ~ 0.002) at 
the scales accessible to our observations (> 1 "). 

4. We find that there can be no systematic differ- 
ence in the non-thermal linewidths between col- 
lisionally excited and recombination lines. This 
contradicts a recent claim that the recombina- 



tion lines are significantly broader (O'Dell et al 



2003) 



from young stars. 
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APPENDICES 

A. SOME PROPERTIES OF THE VELOCITY 
MOMENTS 

The velocity moments of an emission line, as 
defined in equation (JlJ, are all linear functions of the 
line intensity I(V). This means that when two or 
more optically thin emission regions (A, B, C, . . . ) 
are superimposed along the same line of sight then 
the moments of the combined emission can be found 
by simple addition of the moments of the individual 
regions: 



■M n 



M = M ,a + M , B 

■ M 1>B + M 1>c - 
- M 2 n ■ 



Mi = M hA 
M 2 = M 2 , A 



- M 2 ,B 



(Al) 
(A2) 
(A3) 



The first three formal moments M/. are related to 
the more physical quantities: 

total intensity: / = M , (A4) 
mean velocity: (V) = Mi/M , (A5) 

rms velocity width: a = (M 2 /M - (V) 2 ) 1/2 . 

(A6) 

These too can be calculated for the superposition 
of multiple optically thin regions. In particular, for 
the superposition of two regions, A and B, one can 



combine equations (All-(A6l to find 



(V) = 



I = I A + I B , 
(V) a Ia + (V) b Ib 



T B*B 



Ia + Ib 
IaIb((V) a 



{V)b? 



Ia + Ib 



{I a + Ib) 2 



(A7) 
(A8) 

(A9) 



B. FINE STRUCTURE COMPONENTS OF 
HYDROGEN BALMER LINES 

The fine structure components of the Balmer lines 
have velocity differences of < 7 km s _1 , which are 
too small to cause observable line splitting at typical 
nebula temperatures, where the Doppler broadening 
is of order 20 km s _ . However, they do make a 
significant contribution to the observed line width 
and, furthermore, must be considered in order to 
make an accurate determination of the line-center 



wavele ngth. We use the calculations of Cleg g et al. 
( 1999), for densities n < 10 4 cm -3 and interpolated at 



a temperature of 9100 K. The resulting fine structure 
contribution to the velocity dispersion, af s , and rest 
wavelengths, Ao, are given in the first two columns 
of Table H 

Many previous studies have used outdated and 
incorrect values for the rest wavelengths, with a value 
of 6562.82 A being commonly assumed for Ha, rather 
than the correct value of 6562.7910 A. This leads to 
systematic errors of order 2 km s _1 in the derived 
radial velocities of the Balmer lines, as can be seen, for 
example, in Fig. 6 of iBaldwin et al. (20001. However, 



as is shown in Table |4j these discrepancies disappear 
when the correct rest wavelengths are used. 
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